Economic production quantity models for imperfect product 

and service with rework 



o 

(N 



> 
CO 



Allen H. Tai 

Department of Applied Mathematics, The Hong Kong Polytechnic University, Hung Horn, Hong 

Kong 



Abstract 



Oh 

55; 

When imperfect quality products are produced in a production process, rework may 
00 \ be performed to make them become serviceable. In an inventory system, items may 

deteriorate. Selling deteriorated items to customers will create negative impact on 
corporate image. In this paper, two economic production quantity (EPQ) models 
are proposed for deteriorating items with rework process. A single production-rework 
plant system and a system consists of n production plants and one rework plant are 
considered. Approximated analytic results are obtained and numerical examples are 
provided to illustrate the solution procedure. 
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1. Introduction 

Serving high quality products and providing good service can always attract cus- 
^ \ tomers and keep them coming back. However, in reality, production processes are 

often imperfect. For economic and environmental reasons, imperfect quality items 
are reworked to become serviceable again. Due to unsuitable inventory condition or 
other reasons, items stored in inventory face deterioration. In order to provide good 
service, inspection may be carried out to screen out deteriorated items. However, 
when inventory level is huge, full-scale inventory inspection is not possible and only 
part of deteriorated items can be screen out. The remaining deteriorated items will 
then be sold to customers. This will lower customer satisfaction and is harmful to 
corporate image. 

In this paper, we consider two models in which the products produced are imperfect 
and the service provided to customers is also imperfect. The products are imperfect 
in two ways. Firstly, during the production process, imperfect quality items may 
be produced. Good quality items are stocked and sold to customers immediately. 
Imperfect quality items are stocked separately and scheduled for rework. Secondly, 
items may be deteriorated in the inventory. Hence, inspection processes are carried 
out in the inventory. The service is imperfect in two ways. Firstly, inspections for 
deteriorated items in inventory are imperfect. Deteriorated items may be sold to 
customers. Secondly, shortage is allowed and unsatisfied demands are backlogged. 
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Since not all customers accept late delivery, partial backlogging is also considered in 
the paper. 

In the first model, we consider a single production plant system. The plant is also 
capable of rework processing. Part of the imperfect quality items are recovered and 
ready to be sold to customers. At the beginning of the production process, the amount 
of backlog is made up. The production process continues until the total number of 
items produced reaches the economic production quantity Q*. Then rework process 
starts until all imperfect quality items are processed. Demands are then satisfied by 
the serviceable items in the inventory. After the inventory level becomes zero, further 
demands are backlogged and satisfied at the beginning of the next cycle. 

In the second model, we consider a system consists of a central rework plant, which 
is capable of handling imperfect quality items for rework only, and n local production 
plants. The behaviour of the local production plants is the same of the previous 
model expect they can not handle rework process. After the production processes in 
the local plants, the detective items are aggregated and shipped to the central plant 
for rework. The imperfect quality items are then recovered and the recovered items 
are for satisfying the demands at the central plant. At the end of the cycle if there is 
stock at the central plant, they will be sold as one lot at a lower price. The remainder 
of this paper is organized as follows. In Section 2, we give a literature review and the 
motivation of this study. In Section 3, we consider a model for the a single production 
plant. In Section 4, we consider an aggregated model for a central rework plant and n 
local production plants. We then give two numerical examples in Section 5. Finally, 
concluding remarks are given in Section 6 to conclude the paper. 



2. Literature Review 



Economic production quantity (EPQ) is one of the main research topic in pro- 
duction and inventory management. By using EPQ model, optimal quantity of items 
produced can be obtained. Classical EPQ model was developed under various as- 
sumptions. Since then, researchers have extended the model by relaxing one or more 
of its assumptions. 

It was assumed that the items produced is of perfect quali ty in the classical model . 
However, imperfect quality items may be produced in reality. ISalameh fc Jaberl (120001 ) 
proposed an EPQ model with imperfect quality product s. The d efective items are 



screened out and sold as a single batch at a low er price. IWee et al.l ( 120071 ) extended 
the model by considering random defective rate. iJaber et al.l (120081) as s umed the per- 
centage defective per lot reduces according to a learning cu rve. IChand tooi) a pplied 
fuzzy sets theory for modeling defective rate and demands. iRezaei fc Davoodil (120081 ) 
considered a supply chain with multiple products and multiple suppliers. Received 
items from suppli ers were not of per fect quality and decision was made by using 
genetic algorithm. IChung et al.l (120091) proposed an invent ory model with two ware- 
houses, where one of them was rented. lYassine et al.l ( 120121 ) considered disaggregating 
the shipments of imperfect quality items in a single production run and aggregating 
the shipments of imperfect items over multiple pro duction runs. A r eview on EPQ 



models for imperfect quality items can be found in (IKhan et al.l . l201lh . 



2 



How to handle imperfect quality items is another important i ssue. One possible 



ft 

way is to perform rework and make them become serviceable. iChan et al.l (120031 ) 
provided a single EPQ model which considers lower pricing, rework and reject situ- 
ations. The items produced were classified as good items, good items after rework, 
imperfect quality items and rejected items. They assumed that the quality of an 
item was quantifiable and had a normal distribution. IChiu et al.l ( 12004 ) assumed that 
the re work process was im perfect with random scrap rate and shortage was not al- 
lowed. IJamal et al.l ( 12004 ) proposed a model to obtain the optimal batch quality in 
a single-stage production system. Rework was done under two operational policies 
such that total system cost was minimized. IChiu et al.l (120071 ) determined the optimal 
run time for an EPQ mod el with scrap, rework, and stochastic machine breakdowns. 



Buscher fc Lindnerl ( 120071 ) considered a two-stage manufacturing system in which pro- 



duction and rework activities were carried out. The economi c production and r ework 
quantity and the corresponding batch sizes were determined. iTaleizadeh et al.l (1201 if ) 
studied two joint pro duction sys t ems in a form of multiproduct single machine with 



and without rework. lYoo et al.l (120091 ) proposed an EPQ model that incorporated 



both imperfect production quality and two-way imperfect inspection. 

Goods are considered as deteriorating items because their values go down with 
time. Products suc h as electronic produc ts, fashion clothing, food and chemical are 
common examples. iTeng fc Chang] (120051 ) proposed an EPQ model for deteriorating 
items with the demand rat e depended on the selling price of the products and the 
stock level. iLin et al.l ( 120061 ) considered the economic lot scheduling problem (ELSP) 
for deteriorating items. The problem was to schedule multiple products to be manu- 
factured on a single machine repetitively over an infinite planning horizon. iLia ol (l2007h 
developed a production model with finite producti on rate and considere d the effect 
of deterioration and permissible delay in payments. IChung &: Wed (|201lf) considered 



short life-cycle deteriorating items with green product design. IWidyadana fc Wee 
( 120121 ) proposed an EPQ model for deteriorating items with rework, which was pre- 
formed after m production setups. 

One key assumption of classical EPQ model is that no shortage is allow ed. Short- 
age may be handled in two ways: backorders and lost sales. IWee et al.l ( 120061 ) proposed 
an integrated model for deteriorating items in which shortages were completely back- 
ordered. A periodic deliver y policy for a vend or and a production-inventory model for 
a buyer were established. I Wee et al.l ( 120071 ) considered permiss ible shortage backo- 
rdering and the effect of varying backordering cost values. Later, IChang fc Hoi (120101 ) 
used the renewal-reward theorem to derive the exact closed-form solutions of the op- 
timal lot size, backorder ing quantity and maximu m expected net profi t per unit time. 



Cardenas-Barronl (120091 ) extended the models in (IJamal et al.l . 12004) by consid ering 



plann ed backorders. Partial backordering is considered if lost sales is allowed. iMak 



( 119871 ) proposed a opti mal production-inven tory policy for an inventory system with 



partial backordering. iPentico et al.l ( 120091 ) investigated the model in (jMaki . 119871 ) 



and r edeve loped it for the EPQ with partial backordering using simpler expressions. 



Wed (119931 ) proposed an economic produ ction policy for de teriorating items with par- 



tial backordering using iterative method. iGiri et al.l ( 120071 ) considered an EPQ model 
with increasing demand rate and adjustable production rate while shortage are partial 



3 



backlogged. iTeng et al.l (120071 ) gave a comparison between two pricing and lot-sizing 
models with partial backlogging and deteriorated items. 

We notice that not much studies considered a model with imperfect quality and 
deteriorating items, rework and shortage. On the other hand, the effect of selling 
deteriorated items to customers has not been addressed fully also. In this paper, we 
aim at providing analytic results to address to above issues. 



3. The Basic Models 

In this section, we consider a single production plant system which is also capable 
of handling imperfect quality items produced during the production process. The 
following notations are used throughout the paper. 

p production rate (unit/unit time) 

a percentage of good quality items produced 

A demand rate (unit /unit time) 

9 percentage of items deteriorated per unit time 

7 percentage of deteriorated items screened out from the inventory 

p r rework process rate (unit/unit time) 

a r percentage of imperfect quality items recovered 

(3 percentage of customers who accept backlogging 

K setup cost for a cycle ($) 

c deterioration cost ($/unit) 

Cd penalty cost of selling deteriorated items to customers ($/unit) 

c p cost of unrecoverable imperfect quality items ($/unit) 

c s shortage cost ($/unit/unit time) 

c u unsatisfied demands penalty cost ($/unit) 

h s holding cost of serviceable items ($ /unit /unit time) 

h r holding cost of imperfect quality items ($/unit/unit time) 

We have the following assumptions in the development of the model. 

1. The rates p, A, p r and the percentages a, 7, a r , (3 are known constants. 

2. Only serviceable items deteriorate with constant rate 9. 

3. Shortages are allowed and are backlogged. 

4. Backlogged demands are made up at the beginning of the cycle. 

5. Deteriorated items and unrecoverable imperfect quality items are disposed. 

6. Recovered imperfect quality items are considered as good quality items. 

The behaviour of the inventory level of serviceable items at any time during a 
given cycle is illustrated in Figure 1 and the inventory level of imperfect quality items 
at any time during a given cycle is illustrated in Figure 2. 
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Figure 1: Inventory level of serviceable items. 




3.1. Model for Partial Backlogging 

The inventory level of serviceable items at time t over the five periods in a cycle 
are determined by the following differential equations: 

I' 1 (t 1 ) = ap-\, 0<h<Tx, (1) 

r 2 (t 2 )+ 1 ei 2 (t 2 ) = ap-X, 0<t 2 <T 2 , (2) 

I' 3 (t 3 ) + j9h(t 3 ) = a rPr -A, < t 3 < T 3 , (3) 

I'M + jOh(U) = -A, < t 4 < T 4 , (4) 

I' 5 (t 5 ) = -fiX, 0<t 5 < T 5 . (5) 

With the boundary conditions li(0) = J 5 (T 5 ) = -I b , J^Ti) = I 2 {0) = 0, I 2 {T 2 ) = 
J 3 (0) = I s , I 3 (T 3 ) = ^(O) = I m , 74(14) = h(0) = 0, the solutions for the above 
differential equations are 

h(h) = (ap - \)t! - I b , 0<ti<T 1} (6) 

/a(*a)= (^^)(l-exp(- 7 ^ 2 )), < t 2 < T 2 , (7) 

h{h) = (is - exp(- 7 ^ 3 ) + (^^) , < h < T 3 , (8) 

J 4 (t 4 ) = (j m + A) exp(- 7 ^ 4 ) - (A) , < U < T A , (9) 

I 5 (t 5 ) = -fiXt 5 , 0<t 5 < T 5 . (10) 

Hence, it can be deduced from Eq. ([7]) that at t 2 = T 2 , the inventory level of serviceable 
items is 

I s = (^=A)(l-exp(- 7 er 2 )) ) (11) 
and from Eq. that at t 3 = T 3 , the maximum inventory level serviceable items is 

I m = ( QrPr ~ A ) ( X " ex P("7^ 3 )) + exp(- 7 0T 3 )I s , (12) 
and from Eq. ([9]) that at £4 = T4, I^T^) = 0, which gives 

7 -= (A)(exp(70T 4 )-l). (13) 

Also, from Eqs. fl6]) and (TTO]) the unfilled order backlog is 

I b =(ap-X)T 1 = (3XT 5 . (14) 

For the imperfect quality items, the maximum inventory level is given by 

I c = (l-a)p(T 1 + T 2 )=p r T 3 . (15) 

The total operating cost consists of the deterioration cost, penalty cost of selling 
deteriorated items to customers, holding costs of serviceable and imperfect quality 
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items, setup cost for a cycle, cost of unrecoverable imperfect quality items, shortage 
cost and unsatisfied demands penalty cost. Our aim is to minimize the total cost per 
unit time, which can be expressed as 

TC = ~((ap - A)T 2 - J a ) + ~((a rPr - A)T 3 - (J m - Q) + |(J m - AT 4 ) 



+ 



1-7 



7 



|((ap - A)T 2 - I,) + ^((a rPr - A)T 3 - (J m - I.)) + |(J m - AT 4 ) 



+ 



T 



( [ h(t 2 )dt 2 + [ 3 I 3 (t 3 )dt 3 + [ I 4 (U)dt 
v jo Jo Jo 



fer/ (Ti + r 2 + r 3 )p r r 3 

TV 2 

K c v , . m c s (an — A)Ti 2 c s B\T 5 2 c u 
+ - + ^(1 - a r )p r T 3 + ^ * 2 > 1 + + ^/3'AT 5 , 



(16) 



where ft' = 1 — /3. 

In order to simplify the expression, the relationships between Tj (z = 1, . . . , 5) and 
T are needed. Firstly, T = T 4 + T 2 + T 3 + T 4 + T 5 together with Eq. (THj) give 



and 



T-i 



ft\ 



ap — ft'X 
ap — A 



(T-T 2 -T 3 -T 4 



(T-T 2 -T 3 -n). 

ap — ft'X 

Then, with Eqs. ( 1T31) and ( fl5l) . Eq. ( JIB]) can be expressed as 



((ap - A)T 2 + (a r p r - A)T 3 - AT 4 ) 



+ 



/>.s 



ap — A 



( 7 #T 2 + exp(- 7 #T 2 ) - 1) 



h_ _ a r p r - A 
1 7 # ( 7 #) 2 

A 



(1 - exp(- 7 0T 3 )) + 



a r p r — A 



+ 



hoy 



(exp(7#T 4 )-l- 7 #T 4 ) 



+ K\(Pr 2 + (l-a)p-p r )T 3 2 - 



T 



2(1 -a)p 



c s (ap — X)ftX 



c u {ap- X)ft'X 
T ap — p'A 



(17) 
(18) 



(19) 



In order to find the optimal cycle length and the optimal time periods Tj (i = 
1, . . . , 5), we use the following method for getting an approximation of TC for which 
the optimal solution can easily be obtained. Such approxi mation is commonly used in 
modelling systems for det eriorating items, see for example (jWed . ll993l ; IWidyadana fc Wed . 



20121 ; lYang &: Wed . 120061 ). Next, we express T 2 in terms of T 3 , T4 and T. From Eqs. 
f lT5|) and f fT7|) . T 2 can be expressed as 



-T 3 + 



f3\ 



where 



From Eqs. (TTTjh flI2j) and (fI3]h we have 



ap — A ap — A 

(ap - /3'A)p. 



(T 4 -T), 



(20) 



w = /3X + 



a)p 



A(exp(7^T4)-l)-(ap-A)exp(-7^T 3 )(l-exp(-7^T 2 )) = (a r p r -A)(l-exp(-70T 3 )). 

(21) 

In what follows, we use the Taylor series approximation under the assumptions that 
7#T 2 , 7#T 3 and 76T4 are small. We adopt that 



expfxj « 

when x is small. We simplify Eq. ( f2~TT) to 



ar 



A 



a r p r — A 



T 4 + 



— A 
a r p r — A 



T 2 . 



(22) 



The derivation can be found in the Appendix. With Eqs. ( 120]) and (l22j) . the first term 
of TC in Eq. ffl~9]) then becomes 



7c+(l-7)Gi/A#T 4 



T V 2 

Under the assumption that 76T2 is small, I s in Eq. ffTTj) becomes 



(23) 



J s = (ap - A) (t 2 - 



7#T 2 



Similarly, it can be shown that the second term of TC in Eq. ([19]) is 



K 
T 



(ap-\)T 2 2 



+ (ap-\)[T 



lOT 2 



n- 



lOT 3 



| (a r p r - A)T 3 2 | AT 4 S 



(24) 



If we further assume that 1 — 7#T 2 /2 « 1 and 1 — 7#T 3 /2 « 1 then Eq. fl24|) can be 
simplified to 



~T 



(ap - A)T 2 



+ (ap - A)T 2 T 3 + 



(a r p r - A)T 3 2 AT4 



(25) 



Hence the expression of TC can be simplified as 
7 c+(1- 7 )q/A0T 4 2 < 



TC 



ha 

T 



T 

(ap-X)T 2 2 (a r p r -X)T 3 2 XT 2 - 
h (ap - A)T 2 r 3 H 1 — 



K \(p r 2 + (1 - a)p-p r )T 3 2 - 



2(1 -a)p 



(26) 



K c p c s (ap - X)ftX 2 
+ — + —(1 - a r )p r T 3 + —7^773 o7Tt( t ~T 2 -T 3 - T 4 ) 



+ 



T T 

c u (ap — X)ft'X 



T 2(ap - ft'X) 



T ap- ft'X 



f— (r-r 2 -r 3 -T 4 ). 



The optimal values of (T 4 ,T) can be obtained by substituting T 2 and T 3 in Eqs. 
(120]) and (J22J) into Eq. (J25J and solving 8TC/dT 4 = and dTC/dT = simultane- 
ously. Such computation can be done by any mathematical software. 

3.2. Model for Complete Backlog ging 

If we assume that all the unsatisfied demand are backlogged, i.e. (3 = 1 and ft' = 0, 
then the total cost per unit time can be expressed as 



TC 



1C +(l-j)c d /X9T 4 



T 



+ 



h x 



[ap - X)T 2 2 (a r p r - A)T 3 2 AT 4 2 ^ 
h [ap - X)T 2 T 3 H 1 — 



+ 



h r 
T 



(p r 2 + (l-a)p-p r )T 3 ' 
2(1 -a)p 



(27) 



with 



, K . C P ( , s rp c 8 (ap- A) A 
+ T + T (1 - a ^ T3+ T^o^ 



A 



(T — T 2 — T 3 — T 4 ) 2 , 



and 



where 



7i = — (T-T 2 -T 3 -T 4 ), 
ap 

T 5 = ^^(T-T 2 -T 3 -T 4 ) 
ap 



" T 3 + -^(T 4 -T), 



ap — A — A 

ap r 



(28) 
(29) 

(30) 



w = A + 



1 — a 
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For T 3 , we use a simplified version of Eq. (|22|) under the assumption that 1 + 76T4/2 
1: 

rp A ap- X 



a r p r — A a r p r — A 



Together with Eq. f l30]) . we have 
T 3 



-T 



(l-a)A 



a; + a r p r — A ap r + (1 — a)a r p,, 
Finally, we can express TC in terms of T4 and T only: 



-T = T)T. 



(31) 



(32) 



h s 
+ — 
T 



+ T 



(AT 4 - (a r p r - \)r]Tf 



2(ap - A) 



+ A77T4T 



(a r p r - A)?? 2 T 2 AT 4 



(p r 2 + (1 - a)p ■ p r )rj 2 T 2 



2(1 -a)p 



if C. 

— + ^(1 - a r )p r r]T 



+ 



A 



T 2ap(ap — A) 



[((1 - r?)(ap - A) + 77(a r p r - A))T - apT 4 ] 



where 



A = /i 6 



AT + BT A + C^- + ^ + D, 



(a r p r - A) V (a r p r - A)?7 2 ] |"(p r 2 + (1 - a)p ■ p r )rf 



(33) 



+ c s 



2(ap - A) 2 
"A((l - r}){ap - A) + 7?(a r p r - A)) 2 



2(1 -a)p 



2ap(ap — A) 



B = K 



Xr] — 



(a r p r - A)A?7] rA((l - rf)(ap - A) + i](a r p r - A)) 



ap — A 
Atf 



- c. 



— A 



(34) 



C = { 1 c+{l- 1 )c d )[ — )+h s 



A 2 



A 

+ - 



apA 



2(ap - A) 



-2(ap-A) 2. 
L> = Cp(l - a r )p r r]. 

The optimal pair can be obtained by the following theorem. The proof can be 
found in the Appendix. 

Theorem 1. TC(T4,T) is strictly convex. The optimal pair exists if 

B < and AAC > B 2 

and is given by 



(T 4 *, T* 



-B< 



K 



C(AAC -B 2 y V 4AC - B 



CK 



(35) 
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We may also deduce from Eqs. ( TT5|) and (1321) that the optimal production time 
and the economic production quantity are given respectively by 




A 



(36) 



ap + (1 — a)a r p 



and 



q*= p t; 



A 



(37) 



a + (1 — a)ct r 



4. The Aggregated Model 

In this section, we consider a system consists of a central rework plant, which is 
capable of handling imperfect quality items for rework only, and n local production 
plants. After the production processes in the local plants, the imperfect quality items 
are aggregated and shipped to the central plant for rework. We assume that all the 
local production plants are the same as the one considered in Section 3 except rework is 
not preformed there. Therefore, the time between two shipments of imperfect quality 
items is T. It is natural to use T as the cycle length for the central rework plant also. 
Besides the notations and assumptions given in Section 3, we have more as stated in 
the following. 
Notations: 

K c setup cost for a cycle at the central rework plant ($) 

c v penalty cost of selling recovered items at the end of the cycle ($/unit) 

h c holding cost at the central rework plant ($/unit/unit time) 

Assumptions: 

1. Shortages were completely backlogged at the local production plants. 

2. The transportation cost for the imperfect quality items is included in K. 

3. The transportation time for the imperfect quality items is neglected. 

4. Rework is processed in no time and all imperfect quality items are recovered to 
good quality items. 

At each of the local production plant, the behaviour of the inventory level of 
serviceable items at any time during a given cycle is illustrated in Figure 3 and the 
inventory level of imperfect quality items at any time during a given cycle is illustrated 
in Figure 4. 

By using the similar arguments in Section 3, we obtain 




(l-exp(- 7 0T 2 )) 



(exp( 7 flT 4 ) - 1), 



(38) 



I b = (ap - A)Ti = AT 5 , 
I c =(l-a)p(T 1 + T 2 ), 
T = T 1 +T 2 + T 4 + T 5 . 



(39) 
(40) 
(41) 
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Ti + T 2 



time 



Figure 4: Inventory level of imperfect quality items at local production plants. 
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Under the assumption that 7#T 4 is small, we get 



Ti 



— (T-T 2 -T 4 ), 
ap 

ap — A 



-(t-t 2 -t 4 ), 



and under the assumption that 7#T 2 is small, T 2 can be approximated by 



To 



A 7 #T 4 2 



T 4 + 



(42) 
(43) 

(44) 



ap — A V 2 

The optimization problem is to minimize the total cost per unit time, which can be 
expressed as 

u ' x > <apT 4 2 + 1 6\T 4 ^ 
>ap 

h r (l-a)p/XT 7 #AT 4 2 \ 2 



TC = n 



7 c+ (1 - 7 )q / A^T 4 2 \ fe fl / A \ /( 
T V 2 / rUp-A/V 



+ 



T 



V ' 2ap 

| X | c s A/ ap-A \ / y apT 4 + 7gAT 4 2 /2 \2- 
T T2\ «p /\ ap — A / . 



(45) 



+ 



T(T) 
T 



where F(T) is the total cost per cycle of the central rework plant which is defined 
below. 

By using the similar arguments in Section 3 again, the inventory level after rework 
at the central rework plant is given by 



J 6 (t 6 ) = (nl c + A) exp(- 7 W 6 ) -A < t 6 < T 6 , 



A 



7# 
A 



7# 



(46) 



7# 

= n/ c (l - 7 0t 6 ) - At 6 . 
We also have ieC^) = 0, which means 



nT 



(47) 

According to the sale situation, we calculate the total cost per cycle under the 
following two cases, which are illustrated in Figures 5 and 6 respectively: 

• Case I: T 6 > T, i.e. nl c (l - 7#T) - AT > 0. 

In this case, there are recovered items left in the inventory at the end of the 
cycle. The total cost per cycle of the central rework plant is given by 

r-T 



F(T) = h c [J J 6 (t 6 ) dt^j + c v (nl c (l - 7 #T) - AT) + K c 



(4£ 



h, 



nLT - 



+ c v (nI c (l- 1 6T)-\T) + K c . 
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nl. 



time 



Figure 5: Inventory level of recovered items at central rework plant (Case I). 



nl r 




time 



Figure 6: Inventory level of recovered items at central rework plant (Case II). 

• Case II: T 6 < T, i.e. nl c {l - 7#T) - AT < 0. 

In this case, the inventory level of recovered items become zero before the end 
of the cycle. The total cost per cycle of the central rework plant is given by 

F(T) = h c(JJ h{U) d* 6 ) + c u \(T - T 6 ) + K c 

\nI c T 6 ~ ( X + n 2 Icl9 )n 2 ] + cA(T - T 6 ) + K c 



= h r 



(49) 



= K 



+ c u X(T - T 6 ) + K c . 



.2(X + nI c ^6). 

The maximum inventory level of the recovered items is nl c , which can be found 
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by Eqs. (110]), 012]) and 



nl r 



n\(i-a) ( T + i_erl 



a 



(50) 



Now the total cost per unit time can be expressed in terms of T 4 and T only. 
However, the expression is complicated and analysis can not be done easily. Using 
the fact that 7#T 4 and ^OTq are relatively small, we may express Eq. ( I4"5l) in a more 
simple way. Firstly, Eqs. ( 144")) . (T4"T)) and ( 150)) become 



* ^ r ^ nI c T _ nA(l - a) 
-i 4 , ~ -r, "Jc ~ i 

ap — A A a 



and 



n/ c (l - 7 ^T) - AT > i (l - Q 

70 V n(l — a) 



> T. 



(51) 



(52) 



Secondly, Eq. ( |4"5]) becomes 



= n 



2 ^ fc, 
+ 



A 



T Vap — A 



apT 4 



7c+(l- 7 )c tJ / AgT 4 
T V 2 

fe r (1 - q)p /AT\ 2 # c s \/ ap- A \ / 
T 2 lap/ T T2\ ap /I 

F(T) 



apT 4 \ 2 
— A 



(53) 



+ 



T 



Thirdly, F(T) becomes 
F(T) 
h 



nX(l — a) A 
2 



a 



T' 



a 



if J_(i _ a 

1 7 # V n(l-a) 

A(n(l - a)T/a) 2 ] ^ _ n A(l - ^ + 



a 



> T (Case I); 



a / 



7$ V n(l — a) 



< T (Case II). 

(54) 



Since F(T) is a piecewise function, we tackle the optimization problem by considering 
two separate objective functions, which are corresponding to Case I and Case II 



15 



respectively: 



Td = n 



7c + (1 — 7)c d (\6T±\ h s / A \ /apT 4 



T V 2 



+ 7^ 



T Vap — A 



A 



h c 
+ — 
T 



raA(l-a) X Vr2 



a 



ap 

|(=^^(T-T«*)-Ar) + f 



apT 4 
— A 



(55) 



T 4 2 nif + AT C 



AiT + 5T 4 + C-^r + r 



TC 9 = n 



7C+ (1 - 7)c d (\6T±\ / A \(apT A 



V 2 



T\ap-\J\ 2 

+ fe r (1 - a)p / AT\ 2 + X + c, A / ap-A \ /, _ apT 4 \2i 
T 2 Vap/ T T2V ap /\ av — A / 



T 



A(n(l — a)T/a;) 2 ] c u /nA(l — a 



]-|(^^(T-7^)-AT) + § 



(56) 



T 4 2 nK + fT c 



A 2 T + J BT 4 + C^+ .. 



+ D 2 , 



where 



n(l — a)A 2 \ /n(ap — A)A\ /nA(l — a) A 

2a 2 p J + H 2a^ J + H a 2 

2\n „,\2 



/nA^l — a) A\ (nX^6(l — a)\ 

V a 2/ " V a /' 
/ n 2 A(l -a) 2 \ 

V 2a 2 J 5 



2a 2 p 



2ap 



B = —c s nX; 



C = (7c+(l-7)c d )(^)+^( 



nXap \ ( nXap 
+ c s 



2(ap-X)J \2(ap-X) 



D 1 = c v 



nX(l — a) 



a 



Do = c„ A 



' u 



n\(l — a) 



a 



> 
)■ 



(57) 



By Theorem 1, it is not difficult to see that both functions are strictly convex. 
The optimal pairs which minimize TC\ and TC 2 , denoted by (T^ 1 *, T 1 *) and (T 2 *, T 2 *) 
respectively, are 



(Tf X 



-B. 



nK + K c C(nK + K c )\ . = ^ 



C{AAiC-B 2 Y y AAiC-B 2 J 



(58) 



Hence, the following procedure of finding the optimal pair is developed: 
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• Step 1. Solve two optimization problems with objective functions TC\ and TC-i- 

• Step 2. If (1/70)[1 - a/(n(l - a))} < 0, then (T 2 *,T 2 *) is the optimal solution 
to the problem. 

• Step 3. If (1/7$) [1 — a/(n(l — a))} < T 1 *, which means T 1 * does not satisfy 
the condition of Case I, then the optimal pairs must lie on the boundary of the 
feasible region. Hence we replace T 1 * by (1/76*) [1 — a/(n(l — a))} and T 4 * by 
(-B/2C)(lM)[l-a/(n(l-a))}. 

• Step 4. Similarly, if (1/7#)[1 — a/(n(l — a))] > T 2 *, which means T 2 * does not 
satisfy the condition of Case II, then replace T 2 * by (1/76>)[1 — a/{n{l — a))] 
and Tf by (-B/2C)(1/70)[1 - a/(n(l - a))}. 

• Step 5. Compare TCi(T 4 *, T 1 *) with TC 2 (T 2 *, T 2 *). If TC^T}* ,T 1 *) is greater 
than TC*2(T|*, T 2 *) then the optimal solution to the problem is (T 4 *,T*) = 
(Tl\ T 1 *). Otherwise, the optimal solution to the problem is (T 4 *, T*) = (T 2 *, T 2 *) 

We may also deduce from Eqs. ( )42|) and OH]) that the optimal production time and 
the economic production quantity of the local production plants are given respectively 
by 

t; = t* + t; = ^- (t* + ^y~) (59) 

and 2 

Q* = P T; = ^(r* + ?^y (60) 

5. Numerical Examples 

In this section, we provide numerical examples for the models developed in the 
previous sections. We first give an example for the single production plant model in 
Section 3. We assume that p = 6000, a = 0.7, 9 = 0.1, 7 = 0.6, A = 1000, p r = 4000 
and a r = 0.6. For the operating costs, we assume that c = $40, c p = $30, c s = $200, 
c d = $100, h s = $5, h r = $4 and K = $300. We follow from Eq. that the optimal 
pair is given by 

(T 4 *, T*) = (0.1996,0.2891). 
We find T* and T* by solving Eqs. (f22l) and (|30|): 

T; = 0.0519, T 3 * = 0.0247. 

By Eqs. (128]) and (129"]) . we have 

T* = 0.0031, T 5 * = 0.0098. 

From Eqs. (136]) and (137]) the optimal production time and the economic production 
quantity are given respectively by 

T* = 0.0550 and Q* w 330. 
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From Eqs. (TTTT) and ({TBI the maximum inventory is 

Ira ~ 201 

while the maximum inventory during production period is 

I s m 166. 

From Eqs. (fT4"j) and f lT5|) the amount of shortages per cycle is 

J 6 « 10 

while the maximum inventory level for imperfect quality items is 

I c « 99. 

The corresponding optimal total cost per unit time can be found by using Eq. ( fT9|) : 

TC* = $5837.6. 

We remark that although (T 4 *,T*) = (0.1996,0.2891) is the optimal pair to the 
problem with objective function (|33|) rather than the original problem with objective 
function (fl9|) . it still gives a reasonable approximation to the original optimal pair. In 
fact, using numerical method, the original optimal pair is (T 4 *,T*) = (0.2449,0.3493) 
with TC* = $5800.7, which means the approximation error is less than 1%. 

We next give an example for the aggregated plant model to demonstrate the use 
of the procedure developed in Section 4. We assume the same parameters as the 
previous example with n = 5, K c = $250, c v = $10, h c = $3. We follow from Eq. fl58|) 
that the optimal pairs are given by 

(Tl\ T 1 *) = (0.2469, 0.3337), (T 4 2 *, T 2 *) = (0.2248, 0.3038). 

We next check whether (1/ 7 0)[1 -a/(n(l -a))] < T 1 * and (l/^9)[l-a/(n(l-a))} > 
T 2 *. We have 



which means we need to replace (T 4 2 *,T 2 *) by 

(T 4 2 *,T 2 *) = (6.5760,8.8889). 
Finally we compare TCi(T 4 *, T 1 *) = $21917 with TC 2 (T 2 *, T 2 *) = $145868, therefore 



(T 4 *,T*) = (0.2469,0.3337). 
From Eq. (150 p the maximum inventory level for recovered items is 

nl c w 719. 

From Eqs. ( )59|) and (|60|) optimal production time and the economic production 
quantity of the local production plants are given respectively by 




we choose 



T* = 0.0799 and Q* w 479. 
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6. Concluding Remarks 



In this paper, two economic production quantity (EPQ) models were proposed for 
deteriorating items with rework process. We first developed a model for a single pro- 
duction plant system which is also capable of processing rework. We then developed 
an aggregated model for a system consists of n local production plants and one central 
rework plant. The approximated optimal cycle length and the EPQ were obtained. 
Numerical examples were also provided to illustrate the solution procedure. In this 
paper, the demand rate was assumed to be a constant. An extension to this paper 
can be done by considering time- dependent demand rate or stochastic demand rate. 
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Appendix 

Derivation of Eq. p2]> . We have (Eq. (121)) 

A(exp( T 0T 4 )-l)-(ap-A) exp(- 7 £T 3 )(l-exp(-7#T 2 )) = (a r p r -A)(l-exp(- 7 0T 3 )). 
By using the exponential series approximation, we have 



^ m + _ (ap — a)(i — ,m + <*f) ( 7 m - ^ 

--(a r p r - X)(j9T 3 - 



,2 



2 



--(a r p r - A)^T 3 — 



Under the assumptions that j9T 2 and 76*T 3 are small, we have 
\(t 4 + -— P~) ~ ( a P ~ A)^ 2 = (a r p r - X)Ti 
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which means 



as desired. 



Oi r p r — A 



Ot r Pr — A 



□ 



Ta k 



Proof of Theorem 1. Differentiating TC(T 4 , T) = AT + BT 4 + C-£- + — + D with 
respect to T 4 and T, we have 

dTC 2CT 4 

dTC t K CTa 2 
= A 



Solving dTC jdT 4 = and dTC /dT = simultaneously, we get 



{TIT*) 



-B, 



K 



,2 



CK 



C(AAC - B 2 ) ' V 4AC - 5 2 

as desired. Since we must have T 4 * > and T* > 0, therefore 

B<0 and 4AOS 2 . 

The Hessian matrix is given by 

f 2C 2CT 4 \ 



if 



r 



T 

2CT 4 2K" , 2CT 2 



2^3 



+ 



T 3 / 



It is easy to see that since 2C /T > and det(if) = 4CK/T 4 > 0, H is positive 
definite and hence TC(T 4 ,T) is strictly convex. □ 
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